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Abstract 

We consider a family of cylindrical spacetimes endowed with angular momentum 
that are solutions to the vacuum Einstein equations outside the symmetry axis. 
This family was recently obtained by performing a complete gauge fixing adapted 
to cylindrical symmetry. In the present work, we find boundary conditions that 
ensure that the metric arising from this gauge fixing is well defined and that the 
resulting reduced system has a consistent Hamiltonian dynamics. These boundary 
conditions must be imposed both on the symmetry axis and in the region far from the 
axis at spacelike infinity. Employing such conditions, we determine the asymptotic 
behaviour of the metric close to and far from the axis. In each of these regions, 
the approximate metric describes a conical geometry with a time dislocation. In 
particular, around the symmetry axis the effect of the singularity consists in inducing 
a constant deficit angle and a timelike helical structure. Based on these results and 
on the fact that the degrees of freedom in our family of metrics coincide with those of 
cylindrical vacuum gravity, we argue that the analysed set of spacetimes represent 
cylindrical gravitational waves surrounding a spinning cosmic string. For any of 
these spacetimes, a prediction of our analysis is that the wave content increases the 
deficit angle at spatial infinity with respect to that detected around the axis. 

PACS numbers: 04.20.Ha, 04.20.Ex, 04.20.Fy, 04.30.Nk 

1 Introduction 

The study of cylindrically symmetric spacetimes has received a lot of attention in general 
relativity JI], [|. Cylindrical gravity possesses a non-trivial field content, and therefore 
provides a suitable arena where one can test and develop methods (such as solutions gen- 
erating techniques and quantisation procedures [||, d) which are capable of dealing 
with the infinite number of degrees of freedom of the gravitational theory. More impor- 
tantly, cylindrical spacetimes have found application in describing (idealized) situations of 
interest in gravitational physics and astrophysics. One of these applications is the analysis 
of the propagation and interaction of gravitational waves. 

The first family of exact solutions corresponding to time-dependent cylindrical waves 
in a vacuum seems to have been obtained by Beck ||. This family was rediscovered by 
Einstein and Rosen in a systematic analysis of spacetimes which represent cylindrical or 
plane gravitational waves [El. The waves considered by Einstein and Rosen display what 
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has been called 'whole-cylinder symmetry' jj], ||, namely, they are not only cylindrically 
symmetric, but also linearly polarized, so that the metric can be written globally in a 
diagonal form Q. The most general vacuum solution describing cylindrical waves was 
studied independently by Kompaneets and by Ehlers et al. ||. These waves can be 
described by a gravitational model whose configuration space has two field-like degrees 
of freedom which are subject to generalized wave equations [IU|. The general solution to 
such dynamical equations is not explicitly known. 

A different type of physical phenomena that can be associated with cylindrical space- 
times are cosmic strings pT]| . Cosmic strings are topological defects (characterized by a 
non-trivial homotopy group 7i"i) that could be formed during phase transitions predicted 
by grand unified theories [|12[]. In the particular case of a straight cosmic string with a 
non-vanishing mass per unit length, the exterior gravitational field is simply a vacuum 
cylindrical spacetime that presents a conical defect 
sess spin, 



11, 13 



The string may also pos- 

providing the exterior spacetime with an angular momentum |14], The 
existence of straight cosmic strings would have important consequences for astrophysics: 
they would give rise to discontinuities in the microwave background [I(J and act as grav- 
itational lenses [|I7[]. It has also been argued that cosmic strings might have produced the 

Nevertheless, recent observations 
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density fluctuations that led to galaxy formation 
indicate that the anisotropy of the cosmic microwave background has most probably orig- 
inated in an inflationary scenario 



19 



although the possibility that cosmic strings are 
partially responsible for the formation of structure is still open [pOfl . 

In contradistinction to the properties displayed by cylindrical waves, whose spacetime 
is regular everywhere, a straight cosmic string can be described by a cylindrical spacetime 
whose symmetry axis is singular (quasi-regular, to be more precise |21|| ). This singularity 
corresponds to the line source provided by the string [22, |23|. The exterior is a conical 



geometry that possesses a deficit angle and an angular momentum which are proportional, 



respectively, to the mass and spin of the string |L4|, [L5 



Another example of cylindrical spacetimes are those that represent cylindrical waves 
surrounding a straight cosmic string. The analysis of these spacetimes allows a simplified 
discussion of the interaction between gravitational waves and strings. A few (parameter- 
dependent) families of solutions of this type have been found explicitly 0, [25|, |26| . 
Essentially, these solutions have been constructed by applying the soliton technique of 
Belinskh and Zakharov, which makes use of the existence of two commuting Killing vectors 
P, 27j. The spacetimes obtained in this way generally present, at spatial infinity and/or 
around the line source, a dislocation in the direction of the axis [ ]23| , Exceptions are 
the Garriga-Verdaguer solutions [26|, whose metric is diagonal, and a subfamily of the 
spacetimes discussed in [|~~ 



\. A property that is shared by all of these solutions, as well 
as by the gravitational cylindrical waves and the straight cosmic string without spin (but 
not by the spinning string), is that the surface spanned by the rotational and translational 
Killing vectors admits an orthogonal surface, i.e. that the considered isometry group is 
orthogonally transitive [|IJ . In practice, this implies the vanishing of the metric components 
that mix the radial coordinate or the timelike one with the coordinates that describe the 
Killing trajectories. 

One of the authors has recently considered all possible cylindrically symmetric space- 
times that are solutions to the vacuum Einstein equations outside the axis, allowing them 
to contain a non-zero angular momentum p3]. Note that cylindrical waves and straight 
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cosmic strings are included in this family of solutions. For all of the spacetimes in this 
family, it has actually been shown that, once the radial coordinate is chosen orthogonal to 
the group trajectories, the shift vector has non-vanishing projection on the Killing orbits 
provided that the angular momentum differs from zero. Therefore, orthogonal transitiv- 
ity does not hold in this case. In addition, it has also been shown that, in any of these 
solutions, there exists a dynamically conserved quantity that describes the total energy 
per unit length in the axis direction |29j . A similar energy density was known to exist for 
purely gravitational cylindrical waves [TO, BO, FJTJ]. The energy density for those waves is 



actually a non-polynomial function of the total C-energy introduced by Thorne ||, and 
turns out to be positive and bounded from above. Although the total C-energy ceases to 



be a constant of motion when the angular momentum differs from zero, the results of 
generalize the expression of the energy in the presence of an axial singularity endowed 
with spin. In particular, the energy density continues to be bounded both from below 
and above. 

In this paper, we will argue that, apart from a possible dislocation, the spacetimes 
considered in p9| can, in fact, be interpreted as describing the interaction of a spinning 
cosmic string with a cylindrical gravitational wave. In this sense, it is worth noting 
that the degrees of freedom outside the symmetry axis are just those corresponding to the 
cylindrical reduction of vacuum general relativity, namely, the field-like degrees of freedom 
of cylindrical waves. Furthermore, we will prove that the singular behaviour allowed on 
the axis produces a conical geometry which is endowed with a constant angular momentum 
and whose deficit angle does not vary in time. These are precisely the effects of a spinning 
cosmic string. So the intuitive picture that one gets is that of a stringy defect surrounded 
by an empty cylindrical spacetime which differs from the Minkowskian vacuum in that it 
generally contains time-dependent gravitational fields. 

In order to attain this picture, we will investigate the form of the metric, both at 
spatial infinity and around the symmetry axis. We will show that these asymptotic 
metrics can indeed be understood as being created by a spinning cosmic string interacting 
in a non-linear way with a cylindrical wave. We will also see that the main effect of the 
gravitational content at spatial infinity is to increase the deficit angle caused by the string. 
To reach these conclusions, we will first carefully determine the boundary conditions which 
guarantee that the analysed cylindrical solutions are rigorously defined. These boundary 
conditions are crucial to show that the system possesses non-trivial solutions that are 
physically acceptable and admit a well defined notion of energy density. If such boundary 
conditions did not exist, our whole analysis would simply remain as a formal discussion 
devoid of physical content. Furthermore, it is precisely the knowledge of these boundary 
conditions what will allow us to obtain the asymptotic behaviour of the metric and discuss 
the physical phenomena that can take place in spacetime. 

The paper is organized as follows. In section 2, we summarize the main results of |29| . 
This includes the expression of the metric in terms of the field-like degrees of freedom of 
the system, the dynamical equations that these fields satisfy and the Hamiltonian that 
generates the evolution. In section 3 we discuss boundary conditions that guarantee that 
the metric expressions are well defined. The stability of these conditions is studied in 
section 4. Once a set of dynamically stable conditions has been determined, we check in 
section 5 that they really lead to a consistent Hamiltonian formalism. We then analyse 
the behaviour of the metric at large and short distances from the axis in sections 6 and 
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7, respectively. In section 7, we also discuss the physical interpretation of the solutions. 
Section 8 contains the conclusions and further discussion. Finally, two appendices are 
added. 



2 Spinning spacetimes 

In this section, we briefly review the main results of concerning the Hamiltonian for- 
malism for cylindrical spacetimes in vacuo that, in principle, do not include the axis of 
symmetry. The starting point of this analysis is the Hamiltonian formulation of vacuum 
general relativity for spacetimes with two commuting Killing vectors, one of them trans- 
lational and the other one rotational. It is possible to introduce coordinates x a = {z, 8} 
(with z G M, 9 G S 1 and a = 1, 2) adapted to these Killing isometries so that the metric 
is independent of x a . In addition, we will call x 3 = r > the radial coordinate, denote 
spatial indices with Latin letters from the middle of the alphabet and adopt units such 
that 4G — c — 1, where G is the Newton constant per unit length in the z direction. 

All of the gravitational constraints can be removed from the system by imposing 
a gauge fixing that makes use of the cylindrical symmetry. Let us first consider the 
momentum constraints associated with the Killing vectors, which have the form 7i a = 
—2d r (h ai lV r ) . Here, hij is the induced 3-metric and n u its canonically conjugate momen- 
tum ||29|| . These constraints can be eliminated by requiring that the metric components 
h ar vanish. In this gauge, the solution to the constraints TL a = is given by n ar = /i ab q,/4 
(a, b = 1, 2), where the c a s are two constants related to global properties of the spacetime. 
In particular, if the spacetime is regular everywhere, including the axis, the constants c a 
must vanish. The stability of the gauge conditions, on the other hand, determines the 
value of the components N a of the shift vector. After this partial gauge fixing, the system 
has two constraints and the gravitational degrees of freedom reside in the components h a b 
and h rr of the induced metric. 

In fact, the constants c z and eg describe, respectively, the linear momentum in the z 
direction and the angular momentum contained in the spacetime, both of them expressed 
as linear densities. In order to prove this statement, it suffices to remember that, mod- 
ulo surface terms, the generators of the asymptotic rotations (far from the axis) and the 
asymptotic translations in the z direction both have the integral form / drN l 'Hi 1 where the 
shift N l vanishes in the region r C 1 and tends to 5 l a when r — > oo, with a = 1 for trans- 
lations and a = 2 for rotations |50[. Taking into account the expression of 7i a , a simple 
integration by parts shows that the considered generators are differentiable on the grav- 
itational phase space, and hence well defined, provided that the neglected surface terms 
coincide with the limit of 2h ai W r when r tends to infinity. With our gauge fixing, this 
limit is half the constant c a . Therefore, on solutions to the gravitational constraints, the 
values of the analysed generators, which are the linear and angular momentum densities 
under consideration [|3(], |32|], are given in our system of units by the constant quantities 
c z /2 and c$/2, respectively. In the following, we will allow the possibility that the cylin- 
drical spacetime contains a non-zero angular momentum, but will restrict our attention 
to the case in which c z vanishes. Otherwise, our gauge-fixing procedure would lead to 



metric expressions plagued with divergences that would render them meaningless [29 



The gauge freedom corresponding to the only momentum constraint that remains on 
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the system can be removed by choosing as a radial coordinate the square root of the 
determinant of the metric on Killing orbits, since this metric degenerates on the axis and 
its determinant is supposed to possess a spacelike gradient. The canonical momentum of 
this determinant can then be determined by solving the radial constraint on the gauge 
section. The requirement of stability of the gauge condition, on the other hand, fixes 
the radial component of the shift, which turns out to be proportional to the momentum 
P w canonically conjugate to the metric variable w — In \Jh zz h rr . One can then make the 
shift component N r equal to zero by imposing that P w vanishes as a gauge condition for 
the Hamiltonian constraint. In this way, one eliminates all the gauge degrees of freedom 
from the system. The expression of the metric function w can be obtained by solving the 
Hamiltonian constraint with P w — 0. In addition, the dynamical stability of the gauge 
(i.e. P w = 0, where the overdot denotes the time derivative) determines the value of the 
lapse function, assuming that it approaches the unity in the limit r — > oo. 

After this complete gauge fixing, the system can be described by two canonically 
conjugate pairs of field-like degrees of freedom, which we will call (v, P v ) and (y, P y ). The 
reduced metric can be written in the form 



where 



7 2w+y 



N 2 dt 2 + dr 2 + e y r 2 (d9 + N e dt) 2 + e~ y dz - vd9 + {N z -vN e )dt , (2.1) 



e 2w 
N = 
H = 

and the shift vector is 



4E[r]r 2 



c 2 e + ADr 2 - 2c 2 e r 2 / r ds s~ 3 (E[s] - 1) 



2 

r 



-2w ; 



E[r] = exp Qf H 



^ + ^fe- 2 y + P; + P 2 r 2 e 2y 



(2.2) 
(2.3) 
(2.4) 



/•OO Ac _ 

N Z = CgFoo / -rVE[ S ] 



(2.5) 



In these equations, Wqq and Eoo are, respectively, the values taken by w and E[r] in 
the limit r —>■ oo. To arrive at these expressions, the shift vector has been chosen to 
vanish when r tends to infinity and the value of y in this limit has been set equal to zero 
by an appropriate scaling of r, z, eg and the fields (v,P v ). Finally, the parameter D, 
which determines the sub-leading terms in the metric function e 2w around the axis, has 
been assumed to be a positive constant. This assumption is necessary if one wants that 
the dynamics of the gauge-fixed system is generated (at least formally) by a Hamiltonian 
P^| . We will show in section 5 that the constancy of D > is, in fact, guaranteed by the 
boundary conditions that the fields must satisfy. 

The dynamical evolution is dictated by the equations 



2F 00 P v re 2y - 2w E[r 



V = 2F C 



P 



y „-2w 



r 



E[r] 



(2.6) 
(2.7) 
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(2i 



P y = F^dr (?gre-*°E[r]\ - ^e'^r] [APye 2 * - {d r v) 2 e^] , (2.9) 

which are generated in the reduced system, via Poisson brackets, by the Hamiltonian 
Hr = 1 — e~ w °° . Employing equation fl2.2|), one explicitly obtains 



We are assuming that the integrals involved in these expressions converge; we will return 
to this point in section 3. An additive constant in the Hamiltonian has been fixed by 
requiring that the energy vanishes for Minkowski spacetime (i.e. when D = = 1, 
c e = 0). 

Note that, since the reduced Hamiltonian must be real to define an acceptable time 
evolution, the functional expression (|2.10|) implies that the value of Hr, which provides 
the energy density per unit length in the axis direction, is bounded from above by unity. 
The value Hr = 1 cannot be reached, because when e~ Woc vanishes the metric is ill- 
defined, according to our equations. On the other hand, taking into account that H is 
a non-negative function on the reduced phase space, one can see that E[r] > 1 for all 
values of r > 0. As a particular consequence, Q is non- negative. It is then easy to find 
a lower bound for the energy density, provided that the positive parameter D is fixed. 
One gets that Hr > (1 — vD). The lower bound is reached on solutions with vanishing 
momenta P v and P y and constant fields v and y, because it is only then that H vanishes 
and E[r] = 1. Remembering that the value of y when r —>■ oo has been set equal to 
zero and assuming that there exists no dislocation in the direction of the axis P^, |2T| . 
we are then left with only one spacetime. This spacetime, which minimizes the energy 
density and can therefore be regarded as a background for the solutions with parameters 
D and eg, is precisely the region without closed timelike curves (CTCs) in the exterior of 



a spinning cosmic string j29 



As we have briefly commented above, all the integrals over r G IR + involved in the 
metric expressions must converge; otherwise, our previous calculations would lead to 
physically unacceptable solutions. Supposing that the fields {v, y, P v , P y } are sufficiently 
smooth as functions of r over the positive axis, the convergence can be ensured by imposing 
suitable boundary conditions on the fields, both at r = and at infinity. We will discuss 
this point in section 3. Of course, such boundary conditions must hold at all times and, 
therefore, be preserved by the dynamical evolution. This issue will be analysed in section 
4. In addition, it is clear from the form of the metric that, for real fields, the spacetime 
is Lorentzian with t being the time coordinate if and only if e 2w is positive for all values 
of r. In fact, remembering that H is a non-negative function on phase space and that 
D > 0, one can check from equation (|2.2j ) that e 2w is strictly positive in r > provided 
that e 2w °° is positive ||29|| . This last requirement amounts to the condition 2D > c^Q. The 
only effect of this inequality when cq ^ is to restrict the admissible initial data for the 
fields {v,y, P v , P y }. The reason is that e 2w °° is a constant of motion, because it does not 
depend explicitly on time (D is a constant) and commutes with the reduced Hamiltonian 
Hr = 1 — e~ w °° under Poisson brackets []29|. So e 2u, °° remains positive in the evolution 
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if it is originally positive. Hence, the inequality 2D > c$Q needs only be imposed at 
a given, initial time. Note also that Q vanishes if so do all the fields. Consequently, 
for every value of cq and D > 0, there exists a region on phase space around the origin 
v = y = P v = P y = where the considered inequality is satisfied. 

On the other hand, accepting that the dynamical equations for our fields are valid not 
only in the region r > 0, but also when one approaches the axis, a consistency condition 
that must be verified is the constancy of the parameter D. Finally, if the considered 
spacetimes admit a reduced Hamiltonian formulation, it is necessary that the Hamiltonian 
be not only finite for physical solutions, but also differentiable on the corresponding phase 
space. All of these extra requirements will be analysed in section 5, where we will check 
that they are satisfied as a result of the boundary conditions imposed on the fields. 

So, in summary, the situation is as follows. Although we have obtained a formal 
expression for the gauge-fixed metric and determined the equations of motion that sat- 
isfy the degrees of freedom, the existence of solutions with a rigorously defined metric is 
possible only if we can find suitable boundary conditions compatible with the dynamics 
and the constancy of D. Moreover, these solutions correspond to physically acceptable 
spacetimes only if the restriction 2D > c^Q is satisfied by the initial data. Finally, the 
reduced model possesses a well defined Hamiltonian formalism and a constant of motion 
that provides a valid notion of energy density only if the reduced Hamiltonian is differ- 
entiable. Hence, the importance of proving that the system admits a set of satisfactory 
boundary conditions and analysing their consequences. 



3 Convergence of the integrals 

In the rest of the paper, we will assume that the fields {v, y, P v , P y } are sufficiently smooth 
as functions of r over the strictly positive axis. Then, any possible divergence in the 
integrals over r that appear in the metric must arise either at infinity or around r = 0. 
We also assume that the inequality 2D > c^fl (necessary for w to be real) is satisfied 
and that the value of y when r — > oo has already been set equal to zero by a scaling of 



coordinates, fields and parameters |29| . However, we will not impose yet the vanishing of 
the limits v and of the field v on the axis and at infinity, respectively, although we 
will suppose that these limits are, at least, finite. So, for the moment, we will allow the 
presence in those regions of a dislocation in the z direction. 

In the rest of the paper, we will employ the following notation [33|. For any constant 



number a, the symbol / = o(r a ) when r — > indicates that / is much smaller than r a 
when one approaches the axis, so that the limit of fr~ a vanishes at r = 0. Likewise, 
the notation / = d(r a ) when r — > oo means that the limit of fr~ a vanishes at infinity. 
On the other hand, the symbol / = o(g) when r — > means that there exists a strictly 
positive number e such that the limit of r~ e f/g is zero on the axis. In addition, the 
symbol / = o(g) when r — > oo implies that the limit of r e f jg at infinity vanishes for 
a certain e > 0. In those occasions in which it is clear from the context whether we 
are studying the behaviour of the solutions close to the axis or, in contrast, at spacelike 
infinity, we will employ the abbreviated notation / = o(r a ) and / = o(g), obviating the 
appearance of the limit r — > or r — ► oo. More importantly, denoting by £ any of the 
fields {v,y, P v , P y }, we will assume that the behaviour at infinity and around r = is 
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smooth enough to guarantee that the condition £ = o(r a ) when r — > f (with f = or oo) 
implies that <9 r £ = o(r a_1 ) and <9 2 £ = o(r a ~ 2 ) at f. Similarly, if £ = o(r a ) when r — > f, we 
assume that d r £ = o(r a ~ 1 ) and <9 2 £ = o(r a_2 ) in the considered limit. 

We will also write / = 0(r a ) if there exist positive constants Mi, M 2 , and M 3 such 
that \fr~ a \ < Mi, ^fr 1 '^ < M 2 , and \d 2 fr 2 ~ a \ < M 3 close to the axis r = 0. We will 
then say that the function / is of the order or smaller than r a as r — > 0. In practice, as 
we have commented, we will only consider solutions whose fields are sufficiently smooth 
functions of r, including the limit in which one reaches the symmetry axis. For this kind 
of solutions, our definition of the symbol / = 0(r a ) amounts to the existence of the limits 
of fr~ a , d r fr 1 ~ a , and d 2 fr 2 ~ a when r — > 0. 

With the assumptions introduced above, it is possible to check that our metric ex- 
pressions are well defined for all values of r > provided that Eqo is finite and that, if eg 
differs from zero, the integral Q converges. The first of these conditions, Eoo < oo, implies 
that H = o(r _1 ) on the axis and at infinity. Actually, a sufficient condition for Eoo to be 
finite is that H = o(r _1 ) in the two considered regions. Although weaker conditions are 
possible, we will restrict our discussion to this quite general case from now on. Taking 
into account that H is the sum of four non-negative factors, it is readily seen that our 
condition can be equivalently expressed in the form 

d r v = o(l), d r y = o(r _1 ), (3.1) 
P v = o(r- 1 ), P y = o(l), (3.2) 



both when r tends to zero and to infinity, 
axis, one then obtains that, around r = 0, 

v — v + o(r), 



For sufficiently smooth fields v and y on the 
y = yo + o(l). (3.3) 



Here, yo is the limit of y at r — 0, which we suppose finite. Similarly, far from the axis, 

v = Voo + o{l), y = o{l). (3.4) 

In the case that eg differs from zero, one still has to impose that Q is finite. The 
convergence of the integral at infinity is already ensured by the finiteness of i^o. The 
convergence at r = 0, on the other hand, can be seen to require that H = o(r) on the axis, 
and is satisfied, for instance, if the stronger condition H = o(r) holds. Similar arguments 
to those presented above show that this last condition implies that, when r — > 0, 

v = v + o(r 2 ), y = y + o (r), P v = o(l), P y = {r). (3.5) 

These conditions substitute for equations (|3.2|) and ( |3.3|) on the axis when the spacetime 
contains a non- vanishing angular momentum. 

Let us finally discuss the behaviour of the metric components (with our choice of 
coordinates) in the limit r — > 0. Although the axis is in principle excluded from our 
spacetime, the possible singularities at r = are indeed quite weak, allowing the metric 
to be well defined there. In the absence of spin, eg = 0, it is straightforward to check that 
the metric components are finite on the axis when conditions ( |3.2j ) and ( |3.3| ) are satisfied. 
For spinning solutions, on the other hand, the condition v — Vq + o(r 2 ) when r is small, 
together with the convergence of Q, can be seen to guarantee that iV 2 — vN 6 has a finite 
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limit when r — > 0. Then, the only metric component that can diverge on the axis is the 
diagonal t component. The potentially divergent terms in that component come from the 
contribution of — e y [e 2w N 2 — r 2 (N e ) 2 ]. Given the expression of w and that H = o(r), it 
turns out |29j that this expression has nevertheless a finite limit when r — > 0. Hence, the 
inclusion of spin does not destroy the finiteness of the metric at r = 0. 

4 Dynamical Stability 

Since the metric expressions must be rigorously defined at all instants of time, the bound- 
ary conditions on the fields must be satisfied at every single moment and, therefore, be 
compatible with the evolution dictated by equations ( ^6|) -( ^9|) . Otherwise, the system 
would not admit solutions that respect the conditions imposed at r = and at infinity. 
In particular, this would imply that the metric expressions diverge on the solutions of 
the model, so that they would not lead to acceptable spacetimes. In this section, we will 
analyse whether the conditions introduced in section 3 are dynamically stable and, if the 
answer is in the negative, replace them with stronger conditions that are preserved in 
time. 

Let us first study the asymptotic region far from the axis, r — > oo. Remembering that 
e 2w °° and are positive and finite, it is not difficult to check that equations (|3.2|) and 
( p.4[) , together with our equations of motion, imply that 

v = o(l), y = o(r- 1 ), P v = d(r~ 3 ), P y = ( r - 1 ). (4.1) 

Then, our boundary conditions at infinity are stable, because the time derivatives provide 
subdominant contributions, compared to the leading terms in the fields. In particular, the 
behaviour v = o(l) ensures that is time independent. From now on, we will employ 
the notation to remember the fact that this limit is constant. 

The analysis of the stability around the symmetry axis r = is much more involved. 
The solutions with and without spin must be studied separately, because the factor e~ 2w 
that appears in the equations of motion behaves in a different way: it has a finite limit 
when r — > if eg vanishes, but diverges like c^/(4r 2 ) otherwise. In subsections 4.1 and 4.2 
we will study, respectively, the families of spacetimes with zero and non-vanishing angular 
momentum. 

Before doing this, nevertheless, let us point out that there exists at least an infinite 
family of solutions whose behaviour around the axis is the same in the presence or absence 
of spin. These solutions satisfy much stronger boundary conditions at r = than those 
proposed in section 3. Namely, the fields v — v , y — y , P v and P y , as well as their 
derivatives of any order with respect to r, decrease at r = faster than any homogeneous 
polynomial of r. In other words, calling {ip} = {v — t> , y — yo, P v , Py}, we have that the 
limit of r~ a d™ip vanish when r — > for all non-negative values of the integer numbers a 
and m. Examining the equations of motion ( |2.6| )-( |2~P| ), one can see that these boundary 
conditions are indeed preserved by the evolution. The reason is that the right-hand side 
of those equations are given by terms in which we always find one of our fields, or one 
of their derivatives, multiplied by a factor that diverges, at most, like a negative power 
of r. Therefore, all time derivatives turn out to vanish faster than any positive power 
of r on the axis. This proves the stability and implies that the limits of v and y at 
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r = (i.e. vq and yo) must be time independent on those solutions. Note also that the 
above set of spacetimes contains as a particular subfamily the solutions in which one can 
find a neighbourhood r < tq of the axis (with being a positive number) where the 
momenta P v and P y vanish and the fields v and y are constant. The existence of this 
special, infinite family of solutions was already noted in HJ. Finally, let us comment 
that the rapid decrease of the fields at r = in these cylindrical spacetimes indicates 
a trivial interaction between the vacuum degrees of freedom and the axial singularity. 
Consequently, had we restricted our analysis to just this type of solutions, the physical 
interest of our discussion would be severely limited. 

4.1 Vanishing angular momentum 

We will now discuss the stability of the boundary conditions on the axis when the angular 
momentum vanishes. In this subsection, we will restrict our considerations to solutions 
that admit an expansion in powers of r around the symmetry axis. Obviously, these do 
not include the solutions with rapid decrease of the fields at r = commented above. 

When eg vanish, the factor e~ 2w E[r] that appears on the right-hand side of all the 
equations of motion reduces to the constant parameter D, greatly simplifying the calcu- 
lations. In addition, accepting the existence of power series for the fields, the conditions 
and (IpD translate into v = v + 0(r 2 ), y = y + O(r), P v = 0(1) and P y = 0(r). 
Employing this behaviour and the dynamical equation for v, one concludes that, in fact, 
P v = 0(r) and the limit of v must be constant. We will denote this constant value 
by Vq. Finally, equation ( |2.9| ) and our conditions restrict the field y to have the form 
U — Vo + 0(r 2 ). The limit y , on the other hand, does not need to be constant, but can 
vary in the evolution. 

In conclusion, the dynamically stable boundary conditions on the axis are 

v = v c + (r 2 ), y = y + O(r 2 ), P v = 0(r), P y = 0(r). (4.2) 

In particular, these conditions apply to cylindrical waves, case in which D = 1. Although 
the boundary conditions for these waves have already been discussed in the literature 
||. 10 1, we have included their analysis for completeness. More importantly, it seems that 



the results of (JlOj for waves with general polarization contain some mistakes. Equations 
( |3T2"|) , ( |3.4|) and ( |4.2|) correct previous proposals and provide the behaviour that must be 
imposed on the fields in the asymptotic regions far and around the axis. For vanishing 
fields v and P v , on the other hand, these conditions can be seen to agree (modulo the 
supplementary assumption of expansions in powers of r _1 at infinity) with those imposed 
by Ashtekar and Pierri for linearly polarized waves. 

A more detailed analysis of the compatibility of the boundary conditions and the 
dynamics shows that the power series of the fields must be of the form 



v = v c + £ V m r 2m+2 , P V =Y. P m r 2m+ \ 

m=0 m=0 

00 00 

y = y + ^ » 2m+2 , P v = £ Q m r 2m+1 . (4.3) 

m=0 m=0 

The coefficients of these series (except Vq) are functions of the rescaled time r = F^t 

(where we have used that is a constant of motion in the absence of spin). In fact, 
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let us suppose that we know, at all instants of r, the value of yo and the coefficients 
{V m ,Y m , P m , Q m } for all m < n, n being an integer. It is not difficult to check that the 
power expansion of the equations of motion provides then all the information needed to 
determine the next-order approximation to the fields, i.e. the coefficients with subindex 
equal to n + 1. Moreover, the power series turn out to be determined just by the functions 
of time y and Vq and the constant Vq. Actually, once these coefficients are known, the 
lowest-order contributions in the dynamical equations for v and y allow one to obtain the 
values of Po and Qo, respectively. In addition, with equation ( |2.9| ) and the knowledge 
of yo, Vo and Qo, one can fix Yo- The iterative process outlined above leads then to the 
determination of all other coefficients. 

Finally, let us note that, instead of the functions of the rescaled time yo and Vo, one 
can choose as degrees of freedom in the power series the values of yo and of all the sets of 
coefficients {V m , Y m , P m , Q m } at a fixed, initial time r . These values, together with v^, 
completely determine the fields at To. Given such initial data, the equations of motion, 
which are first-order differential equations, fix the evolution of the fields. 



4.2 Case with spin 

Let us now search for stable boundary conditions at r = when c# does not vanish. Like 
in the previous subsection, we will only consider solutions whose fields {v , y, P v , P y } admit 
a power expansion in r around the symmetry axis. Again, this assumption does not hold 
in the family of spacetimes with a rapid decrease of the fields discussed at the start of 
section 4. In addition, we will restrict our discussion to solutions in which the field v has a 
fixed constant limit at r = 0, which we will call Vq. In contradistinction with the situation 
found for vanishing angular momentum, where the stability of the boundary conditions 
requires that v should be constant, a heuristic analysis of the equations of motion seems 
to indicate that now v may actually differ from zero. However, we will concentrate our 
attention exclusively on solutions whose spacelike helical structure in the vicinity of the 



axis corresponds, at most, to a constant dislocation PSL This includes, in particular, the 



case of a spinning cosmic string [14], pLo[| , in which the dislocation is absent. The possibility 
of vo being allowed to vary in time will be discussed elsewhere. 

With our hypotheses, conditions ( |3.5| ) for the convergence of the integrals appearing 
in the metric expressions become 

„ = ^ + 0(r 3 ), y = y + O(r 2 ), P v = 0(r), P y = 0(r 2 ). (4.4) 



On the other hand, according to equation ( |2.2|) , the factor e~ 2w E[r] that is present in all 
of the dynamical equations diverges at r = like the inverse square of r when eg ^ 0. 
Therefore, the condition that P v should be at most of order r is compatible with the 
equation of motion for this momentum only if v = v q + 0(r 4 ). With Vq being a constant, 
the equation for v then requires that P v must satisfy the stronger requirement P v = 0(r 5 ). 
Concerning the conditions on y and P y , we have, in fact, two possibilities. For a generally 
time-dependent value of yo, equation (|2.7| ) imposes that P y = 0(r 3 ). Then, it turns 
out that the subdominant correction to yo must be at least of order r 6 . The reason is 
that, owing to the time independence of the parameter D, the system must satisfy the 
consistency condition P v d r v + P y d r y = o(r 4 ), as we will see in section 5. This requirement, 
together with the stability of the condition on P y , leads to the result y = yo + 0(r 6 ). 
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The second possibility is that the value of yo is a fixed constant, y%. One would thus 
have y = Vq + 0(r 2 ). In this case the equation of motion for y imposes the condition 
P y = 0(r"). 

We may express these two sets of stable boundary conditions in the symbolic form 
w = ^ + 0(r 4 ), P v = 0(r 5 ), y = y^ + 0(r 4+2K ), P y = 0{r A ~ K ), (4.5) 

where the parameter k can be equal to either plus or minus unity, and y^ generally 
depends on time, whereas y$ is an alternative notation for y$. Note that the possibility 



K — — 1 was not included in the discussions of |2{J. 

With the above boundary conditions, one can see that, on solutions to the equations 
of motion, the power expansions of the fields must be of the following type 

oo oo 

v = 4 + £ V m r 2m+ \ = E Pmr 2m+5 , 

oo oo 

V = y ( o K) + E y m r 2m+i+2K , P y = E Q m r 2m+i ~ K . (4.6) 

m=0 m=0 

Except Vq and y$ X ' = y$, which are constants, all the coefficients in these series depend, 
in principle, on the modified time coordinate r = J F^fjdt. This redefinition of time 
absorbs the common factor that appears in the dynamical equations. 

Actually, all the coefficients in these series can be found if one knows, at all instants of 

(k) 

t, the values of Vq, Vb, and ?/q , as well as the value of Y when k = — 1. The proof of this 
statement is sketched in appendix A. The main difference with respect to the situation 
described when the angular momentum vanishes is that now the coefficients of the power 
expansion are not all functionally independent; instead, they satisfy relations which do 
not involve time derivatives. This issue is also discussed in appendix A. As a result of such 
functional relations, one can, in fact, prove by induction that, at a generic initial time To, 
the series ( [4.6|) can be completely determined if one knows at that moment just the values 
of Vq, y^ and all the sets of the form {V^ n) Qs n , P^ n , Y 3n }, where n is any non- negative 
integer. Note that, since the dynamics is dictated by first-order differential equations, 
the commented collection of coefficients, evaluated at r , provide then all the information 
needed to fix the power expansion of the fields around r = at all instants of time. 

From these comments, it is also clear that a possible procedure to construct admissible 
solutions is the following. At a certain initial time, determine the initial values of the fields 
in a region around r = by fixing the collection of independent coefficients given above. 
Try then to analytically continue such initial values to the whole positive semiaxis. If 
the continuation is possible and satisfies the boundary conditions ( |3.2| ) and ( |3.4| ), use 
such initial data, together with the equations of motion, to arrive at a physical solution. 
Otherwise, employ as initial data the result of a smooth matching between the initial 
values obtained around r = and any initial fields in the region far from the axis that 
satisfy the boundary conditions at infinity. 



5 Consistency of the formalism 

In the previous section, we have proved that there exist boundary conditions that ensure 
that the metric expressions, which had been obtained by means of formal integrations, are 
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meaningful at all instants of time on the solutions of the system. In order to prove that 
such solutions lead, in fact, to physically acceptable spacetimes with a rigorously defined 
energy density, we want to show now that the introduced boundary conditions guarantee 
also that the dynamics is fully consistent and that the system possesses a well defined 
(reduced) Hamiltonian formalism. Let us remind that, in the asymptotic region r — > oo, 
the boundary conditions are given by equations Q3.2| ) and (|3.4|), where the limit of v is a 



constant (i.e. t>oo = v^). On the other hand, close to the axis r = 0, the behaviour of the 
fields is dictated by equation ( |4.2| ) if the angular momentum vanishes, and by equation 
( |4.5| ) if eg differs from zero. Remember also that, in this last case, the constant k can 
adopt the values ±1. In addition, it is worth noting that, although this behaviour on 
the axis was deduced assuming that the basic fields admit power expansions in r, the 
conditions apply as well to the set of solutions with rapid decrease of the fields at r = 
discussed in section 4. So, for all of the considered solutions, it will suffice to prove the 
consistency of the Hamiltonian dynamics when requirements fl4.2|) or (|4.5|) are satisfied. 

We will first prove that the dynamics is compatible with the time independence of the 
parameter D, assuming that the equations of motion remain valid in the limit r — > 0. 
Note that, if this compatibility could not be reached, one would be forced either to admit 
that the system does not possess physically acceptable solutions or to try and introduce 
rather artificial sources on the symmetry axis that could account for the constancy of D. 

The constant D determines the sub leading term in e 2w around the axis r = in the 
presence of spin, and the leading contribution when eg vanishes. In the case with non-zero 
angular momentum, using that H = o(r) for r — > 0, one can check 



e 2 - = lr 2 -^Wo(r 4 ). (5.1) 

Co Cn 



Similarly, when there is no angular momentum, the condition H = o(r _1 ) ensures that 
e 2w = 1/D + o(l). Therefore, the hypothesis that D is constant is compatible with the 
dynamics if and only if dt(e 2w ) = o(r 4 ) when ce ^ 0, and dt(e 2w ) = o(l) when eg vanishes. 
On the other hand,the equations of motion for our system, together with the expression 
of e 2w , lead to |2pJ 

d t (e 2w ) = 2F 00 E[r] (P v d r v + P y d r y) . (5.2) 

Since E[0] = 1, we then conclude that the term P v d r v + P y d r y must be of the form o(r 4 ) 
around the axis when eg differs from zero, and of the form o(l) in the absence of spin. It 
is then easy to check that the respective boundary conditions (|4~5l ) and (^4.2| ) guarantee 
this requirement. 

In order to present a well defined Hamiltonian dynamics, our family of cylindrical 
spacetimes must not only possess a real and finite reduced Hamiltonian Hr, but this 
Hamiltonian must also be differentiable on phase space. Otherwise, Hr would not generate 



a true canonical transformation and, therefore, a valid time evolution fl30| . Note that, if 
that happened to be the case, the system would be missing an acceptable notion of energy 
density. We have already commented that, once the integrals and Q that appear in 



equation ( |2.10| ) are known to converge, the reality and finiteness of Hr amount just to 



the condition 2D > c$Q. This inequality is preserved in the evolution, and can hence 



be regarded as a mere restriction on the initial values of the fields [29]. Moreover, there 



always exists a non-empty region of phase space where this inequality is satisfied. On the 
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other hand, it is proved in appendix B that our boundary conditions guarantee also that 
the Hamiltonian Hr is, in fact, differentiable. Hence, our conditions, together with the 
restriction 2D > c 2 d Vt on the initial values, ensure that the Hamiltonian dynamics is well 
defined and that there exists a meaningful constant of motion that provides the energy 
density. 



6 Metric at spacelike infinity 

Once we have determined the boundary conditions that must be satisfied by the basic 
fields of the system, we can discuss the behaviour of the metric close to the axis r = 
and at spacelike infinity, r> 1, \t\. The knowledge of this behaviour will be essential to 
determine the geometrical properties of our solutions and reach a well founded physical 
interpretation of the vacuum spacetimes under consideration. Let us first analyse the 
region at spacelike infinity. 

Taking into account that when r — > oo the behaviour of the fields is governed, for 
all values of the spin parameter eg, by conditions (|3.2|) and (|3.4j) , with Vqq = being a 
constant, a trivial calculation shows 



2w ^e 2w °° +o(l), 



E[r] = E^l + o(l)], e 
N e = -%e w ™ + o(r~ 2 ), N z -vN 6 = o(r~ 2 ). 

The cylindrical metric adopts then the asymptotic form 



(6.1) 



<lt -<' r "^d0) 4 r 2 d6 2 + e 2w °°dr 2 



ds 2 = [l + o(l)] 

+ [1 + o(l)] {dz + d(r~ 2 )dt - [i4 + d(l)]dey . 



(6.2) 



Disregarding the possible existence of a constant dislocation in the z direction (i.e. taking 
v oo = 0) an d neglecting contributions of the form o(l) in the metric, as well as terms of the 
type o(r 2 ) and o(l) in the diagonal 9 component and the z9 component, respectively, we 
obtain precisely the metric that would be created by a spinning cosmic string. Moreover, 
the density of angular momentum of the studied cylindrical solution, cg/2, coincides with 
the spin (per unit length) of the string that would produce the approximate metric. In 
addition, the reduced Hamiltonian of our spacetime equals the energy density of the string, 
namely, 1 — e~ w °°. Obviously, the deficit angle at spacelike infinity is simply this energy 
density multiplied by a factor of 2n. 

The resemblance in the analysed asymptotic region to the exterior metric of a spinning 
string is considerably enhanced when one restricts one's attention to solutions with 



or 



y 



or 



or 



P, 



or 



(6.3) 



Note that this behaviour guarantees that conditions (|3.2j) and ( |3.4| ) are satisfied. A 
particular subfamily of spacetimes in which these asymptotic restrictions hold is that 
formed by all solutions with fields of compact support, namely, (smooth) solutions whose 
fields v, y, P v and P y vanish in the region r > r 1; with r\ being a positive number. 
From the equations of motion of our system, it is not difficult to check that the compact 
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support of these fields is respected in the evolution, so that the noted kind of solutions 
exist. Employing conditions ( |6.3| ), a detailed calculation leads to the following metric at 
spacelike infinity: 



ds 2 



[1 + o{p- 2 )] 



dt 



dO) + pW + e 2w °°dp z + dz 



+o(p 2 )dzd6 + o(p 4 )dtdz. 



Here p 2 



(6.4) 



_l_ e 2w °°Cg/4: is a new radial coordinate. So we see that the approximation 



in the asymptotic region to the exterior metric of a (possibly) spinning string has been 
improved, with respect to the general case, to the level o(l) in the diagonal 9 component 



and at least to the level o(p 



or 



in the rest of the metric components. 



7 Metric near the axis 

7.1 Vanishing angular momentum 

Since the boundary conditions at r = are different for solutions with and without angular 
momentum, the corresponding behaviour of the metric must be studied separately. Let 
us first analyse the case in which the axis is not endowed with spin. The conditions (|4.2|) 
on the fields imply that E[r] = 1 + 0(r 2 ), and it is then straightforward to see that 



ds 2 = [l + 0(r 2 )]e 



'/ii 



-dt 2 + r 2 d6 2 + ^- + e- 2 y°dz 2 



+ 0(r 2 )dzd9 } (7.1) 



where z = z — Vq9 and t = t/yE^ is a rescaled time. In defining this coordinate, we 
have employed that E^ = De 2w °° is a constant of motion when eg vanishes. Neglecting 
contributions or order r 4 to the diagonal 9 component and of order r 2 to the rest of 
the metric elements, we see that, near the axis, the metric describes the exterior of a 
straight cosmic string, except for the possible periodic structure in the z coordinate and 
the allowed time dependence of the factor e yo , which remains constant for a static string. 

Several comments are in order at this point. First, let us remark that the time coordi- 
nates used in our approximations near the axis and at spacelike infinity differ by a scaling 
that depends on the considered solution. It is not difficult to see that this scaling reduces 
to the identity only for the flat solution with vanishing fields v — Vq, y, P y and P v . This 
solution is possible only if Vq = and y = 0. 

On the other hand, the periodicity of 9 leads to the identification of points (9, z) and 
(9 + 2tt, z — 2ttVq). To avoid this spacelike helical structure, we will assume in the following 
that Vq vanishes. Let us also note that there exist solutions in which y Q is actually constant. 
For instance, this occurs in the solutions with fields that decrease rapidly at r = 0, as 
discussed in section 4. Furthermore, we proved in subsection 4.1 that the time-dependent 
coefficients yo and Vo could be considered as the true degrees of freedom contained in 
the expansions (|4.3| ) around the axis. Therefore, the possibility that y is set equal to a 
constant is indeed available. 

In addition, it is known that (at least when the angular momentum vanishes) cylindri- 
cal gravity in a vacuum can be reduced to three dimensions by employing the existence 
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of a translational Killing vector [[j], 30]. This reduction leads to three-dimensional gravity 



with axial symmetry coupled to two scalar fields, namely v and y [[54]. From this point 
of view, the metric in three dimensions associated with our spacetimes in the vicinity of 
the axis is ds 2 = —dt 2 + r 2 d6 2 + dr 2 /D, modulo relative corrections of the order of r 2 or 
smaller. This is exactly the 3-metric created by a point particle of mass m = 1 — \f~D 
T4J . With this perspective, the time variation displayed by y in equation ( [7.1[ ) can be 



understood just as a time-dependent lift of the metric from three to four dimensions. This 
time dependence is generally necessary to compensate the radial variation of the fields 
around the axis, so as to finally arrive at a vacuum solution for r > 0. 

We also note that the metric ( f^TTf ) presents a deficit angle around the axis, a property 
that is characteristic of stringy defects. Moreover, this deficit angle is independent of 
time, namely 2n(l — VD), as it would correspond to a static string with linear energy 
density equal to m, i.e. the mass of the point particle obtained in the three-dimensional 
version of the system. Thus, we see that, like in the case of a straight string, the effect 
of the singularity on the axis is to introduce a constant deficit when D ^ 1. Note that, 
in order for this deficit to be positive, one must restrict the parameter D to be smaller 
than unity. For D = 1, on the other hand, the spacetime is completely regular E3] and 



the considered solutions are purely gravitational waves. Finally, it is worth commenting 
that the approximations of the metric given in equations ( |6.2|) and (|7. 1|) can be checked to 



describe the behaviour of the solutions with quasi- regular axis analysed in Refs. [24, E 



as well as those contained in the soliton spacetimes of class A2 and B2 constructed by 
Garriga and Verdaguer ||26|| . In agreement with the above discussion, such solutions have 
been interpreted as representing the interaction of cylindrical waves with a cosmic string. 



7.2 Case with spin 

When ce ^ 0, there exist two acceptable sets of boundary conditions (assuming a constant 
value for v at r = 0), each of them leading to a different metric behaviour. We will first 
consider the case in which the field y may have a time- dependent limit on the axis. The 
corresponding boundary conditions are given by equation ( |4.5| ) with k — 1. Remember 
that these conditions apply as well to solutions with fields that decrease rapidly at r = 0, 
even though they were originally deduced for fields that admit an expansion in powers of 
r. One can then easily see that 

E[r] = 1 + 0(r 6 ), 

N e = ^Foo[l + 2fir 2 + 0(r% N z - vN e = N* + 0(r 2 ). (7.2) 
Here, Q is the convergent integral defined in equation ( |2. 1C| ) and 

NS = c e F 00 r^(v-v c )E[r], A = * (7.3) 
Jo r 6 2V-D 

Note that Nq is finite, given our boundary conditions. Introducing then the coordinates 

t = \[D \ l dtF^t), e = e + c e t dmfyF^t), 

Jo Jo 

p 2 = r 2 + A 2 , z = z-v c 6+ [ dtNZ(t), (7.4) 

Jo 



-2w 



4 r 2 



Cylindrical waves and spinning strings 



17 



the metric around the axis r = can be expressed in the form 



ds 2 = [1 + O(r 6 )}e y o } 



(di-Adef + p 2 d9 2 + ^- + e- 2 ^dz 2 



+0(r 2 )dtdz + 0(r i )dzd9 + 0(r 4 )(dt} 2 . (7.5) 

In the definition of our new coordinates, we have taken into account that the quantities 
Fqo, Nq and Q may implicitly depend on time when eg ^ 0. Like in the absence of angular 
momentum, the time and angular coordinates of the expansions around the axis and at 
spacelike infinity [see equation ( |6.2|) 1 differ for all but the flat solution, which has vanishing 
momenta P v and P y and constant fields v = v q and y = 0. Concerning the axial coordinate 
z, we see that the periodicity of 9 introduces again non-trivial identifications of points in 
the sections of constant time unless Vq vanishes. We will thus restrict our discussion to 
the case Vq — from now on. The coordinates z and z can then be seen to coincide on 
the flat solution; otherwise they generally differ. Finally, note that the new coordinate p 
is defined only over the semiaxis (\A\, oo), because, in our spacetime, r must be positive. 
In principle, however, one could try and make an analytic extension to values of p smaller 
than \ A\. This would correspond to imaginary values of r. It is important to remark that, 
in any case, such an extension would give rise to the appearance of CTCs because, from 
equation ( |7.5| ) and at least in the region with p smaller but close to \A\, the diagonal 9 



component of the metric would then be negative, indicating the presence of a timelike 
vector field with closed orbits. 

In terms of p, the expression 0(r 2a ) can be rewritten in the form 0([p — \A\] a ) for 
any real number a. On the other hand, disregarding the small iz and z9 components, the 
contributions of the type 0(r 4 ) to the diagonal t component, and all relative corrections to 
the metric of the order of r 6 , we arrive at a metric near the section p — \A\ that describes 
the gravitational field created by a spinning string, except for the fact that now yo may 
depend on time. Of course, this time dependence does not show up in the solutions with 
rapidly decreasing fields, since yo is then a constant. In addition, from our discussion in 
subsection 4.2, where we showed that the function yo(t) could be considered as one of the 
degrees of freedom contained in the series (f4.6p, it follows that there exist solutions with 



power series at r = in which yo remains constant in time. 

In the general case in which yo depends on time, a line of reasoning similar to that 
presented for spacetimes with vanishing momenta allows us to regard the approximation 
of the metric around p = \A\ as a time-dependent lift to four dimensions of the three- 



dimensional metric produced by a rotating point particle with mass [[TJj]. In addition, 
it is worth noting that the approximated metric again presents a deficit angle that is 
constant. In this spirit, the arguments given in subsection 7.1 support the interpretation 
of the metric in the analysed region as that associated with a spinning string interacting 
with an environment of gravitational waves. This spinning string is characterized by two 
parameters, namely, the corresponding density of angular momentum cg/2 and the deficit 
angle 27r(l — VD). Actually, this value of the deficit is not modified by the presence of 
spin, since it takes the same expression when eg vanishes. 

To end this section, we will discuss the behaviour of the metric around the axis r = 
when the limit of y is a fixed constant y^, so that the boundary conditions are given by 
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equation (|4.5| ) with k = — 1. These imply 




Foofl + 2Qr 2 + 0{r% N z - vN e = N* + 0(r 2 ). 



E[r\ = 1 + 0(r 4 ) 



(7.6) 



Then, employing equation ( |7.4| ), we can write the metric in the form 




(7.7) 



Here, e y = e y o + 0(r 2 ). Note that, around p = \A\, the metric presents, in fact, a constant 
deficit angle. The corrections to this deficit are at most of the order of r 4 , even if the 
subdominant contributions to y$ can be of the type 0(r 2 ). Like in the case k — 1, we will 
restrict our attention to the possibility that Vq vanishes. Since y$ is a constant, it is then 
clear that the metric in the region < p — \A\ <C 1 can be interpreted as that originated 
by a string with spin equal to Co/2 and a conical deficit given by 27r(l — y/~D). 

8 Discussion and conclusions 

We have considered the most general cylindrical solution to the Einstein equations in 
vacuo that, in principle, does not contain the symmetry axis r = 0. This axis, which may 
be singular, has been allowed to possess spin. However, to arrive at metric expressions 
that are well defined in the whole region r > [P9| , we have supposed that the linear 
momentum in the axis direction vanishes. These spacetimes were analysed in a recent work 
P§|] , where a gauge- fixing procedure that removes all of the gravitational constraints was 
introduced. The resulting reduced model can be described by four fields on phase space, 
{v, y, P v , P y }, whose dynamics is generated, at least formally, by a reduced Hamiltonian 
Hr. The singularity on the axis can be characterized by two parameters: a real constant 
cq that provides (twice) the density of angular momentum and a positive number D that 
determines the behaviour of the purely radial component of the metric around r = 0. The 
latter of these parameters should be a constant in order for the dynamics to be consistent 
and the system to possess a conserved energy density PP| , given by the value of Hr. 

We have first searched for boundary conditions which ensure that the reduced system 
obtained after completing the gauge fixing is well defined. These conditions describe the 
behaviour of the basic fields {v , y, P v , P y } at spacelike infinity, r>l, \t\, and near the axis 
r = 0. More explicitly, the consistency of the system implies the following requirements. 
First, all the metric expressions which have been found by means of formal integrations 
must be meaningful. Second, in order for the reduced Hamiltonian to generate a well de- 
fined dynamics and provide the energy density of the system, Hr must be real, finite and 
different iable on phase space. Third, these requirements must be satisfied at all instants 
of time and, therefore, the boundary conditions must be dynamically stable. Finally, 
assuming that the equation of motion for the metric function w (that appears in the 
diagonal radial component) remains valid in the limit r — >• 0, one has to check that the 
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value of the parameter D is, in fact, preserved in the evolution. The existence of bound- 
ary conditions that satisfy the above requirements is most fundamental; otherwise, the 
reduced dynamics would not be consistent and there would not exist physically accept- 
able (non-trivial) solutions in cylindrical vacuum gravity endowed with a non-vanishing 
angular momentum. 

In our discussion, we have supposed that our basic fields are smooth over the whole 
semi-axis r > for all possible values of the time coordinate t. In addition, denoting by £ 
any of these fields, we have assumed that both the requirements £ = o(r a ) and £ = o(r°) , 
either at r = or at infinity, automatically imply that the derivatives <9 r £ and <9 r 2 £ display 
a behaviour similar to that of the field £, but with the exponent a replaced with a — 1 and 
a — 2, respectively. With this assumption, the boundary conditions at spacelike infinity 
turn out to be given by equations (|3.2|) and Q3.4Q . In the latter of these equations, Voo 
(i.e. the asymptotic limit of v) is a fixed number that reflects the possible existence 
of a constant dislocation in the z direction. To avoid the appearance of an asymptotic 



spacelike helical structure [p8[| , one only has to make = 0. 

The boundary conditions that must be introduced on the axis r = 0, on the other hand, 
differ for solutions with or without angular momentum. When eg vanishes, the conditions 
are given by equation (|4.2|). Again, to prevent the existence of a screw dislocation in the 



axis direction |p8fl , one must demand that v vanishes at r = 0, i.e. Vq = 0. For non-zero 



spin, there exists more than one set of admissible boundary conditions. We have only 
considered in detail the case in which the limit of v on the axis r = is constant, thus 
corresponding to a fixed dislocation. Then, the behaviour of the fields depend on whether 
the limit of y is also fixed and time independent or, in contrast, is allowed to vary in time. 
In both situations, the boundary conditions can be written in the symbolic form ( |4.5p , 
where k = — 1 for the first of the considered possibilities (y^ 1 ^ is constant) and k = 1 
otherwise. 

By means of a careful analysis of the equations of motion and the Hamiltonian, we have 
proved that these boundary conditions are stable and that, together with the requirement 
2D > CgQ (which restricts the admissible initial data for the basic fields), guarantee that 
the metric and Hamiltonian dynamics are rigorously defined. We have also checked that 
the evolution is compatible with the constancy of D. So all consistency conditions are 
satisfied. In addition, employing the boundary conditions, we have been able to calculate 
the approximate form of the metric in the region close to the axis r = and at spacelike 
infinity. The physical picture that arises from this analysis is the following. 

In the region near r = (and assuming that Vq = 0), the metric describes a 'rotating' 
conical geometry, with constant deficit angle and a generally non-vanishing angular mo- 
mentum. This approximate metric corresponds to the lift from three to four dimensions 
of the metric produced by a point particle with constant mass and spin given by 1 — \f~D 
and cg/2. The lift is in general time dependent, to account for the possible variation of 
the fields around r = and yield a solution to the Einstein equations at all points of the 
spacetime. The 4-metric can thus be regarded as that caused, in the analysed region, by 
a spinning cosmic string with linear density of energy and angular momentum determined 
by the mass and spin parameters of the analogous particle in three dimensions. Note that, 
in order for this energy density to be non-negative, one must restrict the parameter D to 
be equal or smaller than unity. It is also worth remarking that the time and cylindrical 
coordinates that are naturally associated with this spinning string do not coincide with 
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those selected in our gauge fixing, which are specially adapted to the asymptotic region 
far from the axis. The relation between both sets of coordinates is given by equation (|7.4j) 
(with Vq = 0). This relation continues to be valid even if eg vanishes [see equation (|6.2|) 1. 

Of particular importance is the change of radial coordinate in the presence of spin. 
This change is precisely that which would remove the region with CTCs from the exterior 
of a spinning cosmic string, mapping the resulting spacetime to the sector r > 0. This 
opens the possibility of analytically continuing our spacetime to the region of imaginary 
values of r (namely, to p < \A\) at the price of introducing CTCs. Note that, in our 
original spacetime, the existence of CTCs is actually precluded. An infinite family of 
solutions in which the continuation to imaginary values of r can be straightforwardly 
performed is that with fields of rapid decrease at r = 0, considered in the beginning of 
section 4. It is clear that the metric of any of these spacetimes can be matched smoothly at 
r = with the metric of a spinning cosmic string describing the region p < \A\, providing 
in this way an extended solution that contains timelike orbits generated by dg. 

Since our spacetimes are the most general solution to the Einstein equations in a 
vacuum for r > 0, the only field content outside the axis r = is that corresponding 
to gravitational waves. Consequently, the studied spacetimes represent an ensemble of 
gravitational waves surrounding a singular axis that is characterized by its spin density 
and by the constant deficit angle detected in its vicinity. These are precisely the effects 
that a spinning string would produce. In this sense, the analysed spacetimes can be 
regarded as describing the most general interaction that is permitted in general relativity 
between cylindrical waves and strings with constant density of energy and spin. 

This interpretation is in accordance with the asymptotic form of the metric (|6.2| ) at 
spatial infinity. It is known that, in the region r ^> 1, a cylindrical wave causes a deficit 
angle that is proportional to the total energy density contained in the wave; namely, the 
deficit is 2n{\ — ||10|| . In addition, the wave does not carry angular momentum. 

Therefore, in the case that the gravitational wave surrounds a cosmic string, the deficit 
angle at spatial infinity should be produced by the combined effect of both phenomena, 
and the angular momentum should be originated in the string. So, for a cylindrical wave 
interacting with a spinning string, one would expect the asymptotic metric at spacelike 
infinity to describe a conical geometry with a spin parameter equal to that of the string, 
but with a different constant deficit angle. This is, in fact, the result that we have 
obtained. 

Let us analyse in more detail the relation between the deficit angles encountered at 
r = and at infinity: 2n(l — \^D) and 2tt(1 — e~ w °°), respectively. In the absence of spin, 
the relation is relatively simple; at spacelike infinity, one only has to divide the parameter 
D by -Eqo, a constant of motion that determines the energy of the wave. This quotient can 
be interpreted in terms of the associated C-energy M as an addition of energies. In the 
spinning case, on the other hand, the relation is much more complicated. The quantity 
is not a constant of motion anymore, and the gravitational wave does not possess a 
preserved energy density by its own (since the C-energy is not conserved). Moreover, the 
existence of angular momentum introduces corrections to the deficit angle that depend 
on the gravitational fields not just through the functional E^, but also via Q. The spin 
decreases the effective value of D by an amount of CqQ/2. The resulting deficit is, apart 
from the usual factor of 2ir, equal to the total energy density of the system, given in 
equation (|2.10| ) and which is again conserved. Furthermore, as we commented in section 
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2, the positivity of H on phase space implies that Hr > 1 — yD. Consequently, in 
absolutely all of the spacetimes that we have considered, the deficit angle at infinity is 
greater than or equal to that around the axis r = 0. This is a general prediction that, in 
principle, could be verified experimentally, had we the possibility of making measurements 
in cylindrical gravitational systems containing a stationary stringy defect or, alternatively, 
in systems that could mimic the dynamics of the gravitational field in this situation. 

Our result can also be rephrased by saying that, although the energy of the composite 
system is not the sum of the energies corresponding to the spinning string and the cylin- 
drical wave, the presence of a gravitational wave always results in an increase of the total 
energy. It is not difficult to check that the only situation in which the two studied deficit 
angles coincide is when the fields v and y are constant and the momenta P v and P y vanish 
all over the spacetime, i.e. for the flat solution which (modulo a constant dislocation, 
given by the fixed value of v) describes the vacuum region, free of CTCs, in the exterior 
of a string with mass and spin densities equal to 1 — \fT) and cq/2. 

It is worth commenting that, although the degrees of freedom of our family of space- 
times are those corresponding to purely gravitational cylindrical waves (namely, the fields 
v and y and their canonical momenta), the equations of motion that dictate the evolution 
in configuration space (i.e. in terms of v and y) are not truly hyperbolic partial differen- 
tial equations, except in the absence of spin. The reason is the appearance of the factor 
e~ 2w E[r] in all of the dynamical equations ( |2.6| )- (|2.9| ). This factor has a local dependence 
on the fields if and only if cq vanishes, in which case it reduces to the constant D. In other 
words, the existence of angular momentum introduces, via the process of gauge fixing |29|] , 
a high non-locality in the dynamics of cylindrical gravity. 

Finally, in order to confirm the interpretation that we have put forward for our family 
of spacetimes as cylindrical waves interacting with spinning strings, it would be interesting 
to analyse the Riemann tensor of the solutions, paying a particular attention to the 
contributions that, in the form of distributions concentrated on line sources, could account 



for the appearance of the axial singularity |14|, 35|. In addition, one might also study the 



behaviour of the Riemann tensor at null infinity, determining the radiative content of the 
gravitational field in this region (as it was done, e.g., in Refs. [^5], ^6j). These issues will 
be the subject of future research. 
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Appendix A 

In this appendix, we will discuss the freedom available in the choice of coefficients in the 
series (|4.6|) . We will first show that all the coefficients in these series can be determined 

(k) 

from the knowledge, at all instants of time r, of the values of t>Q, Vq, and y^ , as well 
as the value of Y Q when k — — 1. In order to demonstrate this statement, let us call T n 
the set formed by t> , Uo^ an d all the coefficients {V m} Y m , P m , Q m } with m < n, where 
n is an integer. By expanding the dynamical equations (|2.6|) - (^9|) in powers of r, one 
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can check that the coefficients with subindex equal to n + 1 can be determined from the 
knowledge of T n at all values of r. Consequently, all the information needed to fix the 
series is contained in F\. Let us now analyse the freedom available in the choice of this 
set of coefficients. One can check that the lowest-order contributions to d T v and d T y in 
equation ( |2.6|) and (]2.7|) , respectively, determine P and Q in terms of the constant jjq 
and the derivatives of Vq and Y with respect to r when k — — 1, or as functions of t/q 
and the derivatives of this coefficient and Vq if k — 1. In this last case, in addition, the 
contributions of order r 3 to d T P y turn out to fix Y in terms of Vo and , once Qo has 
been found. This concludes the proof of our assertion. 

Let us consider now the collection of coefficients that appear in the series ( |4.6|) at a 
certain initial time To, rather than as functions of time, and discuss the freedom that exists 
in the choice of such initial data. It is not difficult to check that the coefficients of these 
power expansions satisfy functional relations which do not involve time derivatives, so 
that they are not all independent at r . In particular, for k = 1, the lack of contributions 
of order r and r 3 in d T P v implies that V\ and V> are proportional to Vo- I n addition, 
the requirement that the subdominant correction to d T yQ should be of order r 6 turns 
out to fix the coefficients Q\ and Q2 as linear homogeneous functions of Qq. Similarly, 
for k = —1, the vanishing of the derivative of P y with respect to r at orders r and r 3 
leads to a functional dependence of the coefficients Y\ and Y 2 on yfi, Y and Vo. Taking 
into account this dependence and demanding that d T P v does not include contributions 
of order r and r 3 , one can also determine the coefficients V\ and V2 as functions of Y 
and V). In general, more complicated relations appear when higher-order coefficients are 
considered. As a consequence of these relations, one can prove by induction that, in order 
to determine the series (|4.6| ) at r , it suffices to know at that moment the values of Vq, 
l/o^ and all the sets of coefficients { Vs n , Qs n , P^n, Y Sn } with n > 0. 



Appendix B 

In this appendix, we will prove that the reduced Hamiltonian Hr is differentiable on phase 



space once one adopts the boundary conditions (|4.5| ) [or ( |4.2| ) if the angular momentum 
vanishes], (|3.2|) , and (|3.4|) [where Voo = f^J. The variation of the reduced Hamiltonian is 
given by 

5H R = e — dr5H[2 + c 2 E oo F^ ^E[s}\ , (B.l) 




where H is the function on phase space defined in equation ( j2.4Q . The integrals that 
appear in the term in parentheses are all convergent for r > 0, owing to our boundary 
conditions. In addition, from a variation of the canonical momenta P v and P y , one obtains 

5H = ^{P v r 2 e 2y 5P v + P y 5P y ). (B.2) 

Therefore, the Hamiltonian is differentiable with respect to these momenta if and only 
if the integral over r in the expression of 5Hr converges (both at r = and at infinity) 
for all possible variations of P v and P y . Given the asymptotic behaviour ( p.2| ) and the 



conditions on the axis, equation (|4.2|) or equation (|4.5|) , the admissible variations of P v 



and P y turn out, in general, to be of the same order as the momenta themselves. It is 
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then a simply exercise to check that the studied integration over r leads, in fact, to a well 
defined variation of the reduced Hamiltonian. 

The analysis of the variations of v and y is more complicated. In this case, one gets 



5H 



1 r 



AP 2 r 2 e 2y 5y - (d r vfe- 2y 5y + r 2 d r yd r (5y) + d T ve- 2y d r (Si 



(B-3) 



Note that this variation contains the radial derivatives of 5v and 5y. To get rid of these 
derivatives, one must perform an integration by parts. Consequently, the variation of the 
reduced Hamiltonian splits into a surface term and an integral expression. The surface 
term is 



(r 2 d r y5y + d r ve- 2y 5v) h + cjE^ jf -E[s 



4r 



(B.4) 



where we have employed the notation f\ b a = f(b) — f(a). The integral contribution to SH R , 
on the other hand, can be obtained by replacing 5H on the right-hand side of equation 
( |B1D with 



4P 2 re 2y 5y - (d r 



-2y 



-by - d r (rd r y)Sy - d r d r v- 



8v, 



(B.5) 



and adding to the result the factor 



r 2 _ 

dr-^E^FlElr] r 2 d r y5y + d r ve~ 2y 5i 



(B.6) 



The differentiability of the reduced Hamiltonian then requires that the surface term 
( |B.4|) vanishes and that the integrals over r that determine the variation of Hr are conver- 
gent for all possible values of Sv and 8y. Actually, according to the boundary conditions 
( p.4|) (with Voo = v^a), one has that 5v = o(l) and 5y = o(l) when r — > oo. Here, we 
have imposed that the variations of v preserve the value of v^, since this is a given con- 
stant. Similarly, from equation ( |4.2|) one concludes that, when eg vanishes, the acceptable 
variations of our fields display the behaviour 5v = 0(r 2 ) and 5y = 0(1) around the axis. 
Finally, in the presence of spin, one obtains from equation (|4.5|) that 5v = 0(r 4 ) at r = 0, 
whereas the variation of y admits two different types of behaviour, depending on whether 
y^ is a fixed constant or not. In the first case (k = —1), one arrives at Sy = 0(r 2 ); in the 
second case (k = 1), one gets 5y = 0(1), because now it is possible to vary the limit of 
y on the axis. With this information about the variations of the fields and our boundary 
conditions, it is not difficult to check that the reduced Hamiltonian is indeed different iable 
with respect to v and y. We therefore conclude that our boundary conditions guarantee 
that H R is differentiable on phase space, as we wanted to show. 
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